Superfluid toroidal currents in atomic condensates 
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The dynamics of toroidal condensates in the presence of condensate flow and dipole perturbation 
have been investigated. The Bogoliubov spectrum of condensate is calculated for an oblate torus 
using a discrete- variable representation and a spectral method to high accuracy. The transition from 
spheroidal to toroidal geometry of the trap displaces the energy levels into narrow bands. The lowest- 
order acoustic modes are quantized with the dispersion relation lu ~ |m|w s with m = 0, ±1, ±2, .... 
A condensate with toroidal current k splits the \m\ co-rotating and counter-rotating pair by the 
amount: AE w 2\m\h 2 K,{r~ 2 } . Radial dipole excitations are the lowest energy dissipation modes. 
For highly occupied condensates the nonlinearity creates an asymmetric mix of dipole circulation 
and nonlinear shifts in the spectrum of excitations so that the center of mass circulates around the 
axis of symmetry of the trap. We outline an experimental method to study these excitations. 

PACS numbers: 03.75.Kk,03.75.Lm 
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I. INTRODUCTION 



The elementary excitations of trapped Bose Einstein 
condensates have been extensively studied in recent 
years. These collective modes are coherent macroscopic 
matter waves that can be used in many applications in 
cold atom Physics. Since the trap geometry defines the 
mode spectrum and amplitudes, recent studies have con- 
sidered spheroidal condensates IE 12 and in novel topolo- 
gies such as toroidal traps 0, 0, laj ■ A toroidal trap is 
of particular interest as it can be employed as a storage 
ring for coherent atom waves or ultracold molecules 
enabling investigations of persistent currents, Joseph- 
son effects, phase fluctuations and high-precision Sagnac 
or gravitational interferometry. More adventurous possi- 
bilities for toroidal condensates include the construction 
of a mode-locked atom laser y| and the creation of sonic 
black holes in tight ring-shaped condensates A nar- 
row ring of condensate, effectively reducing the dimen- 
sionality to ID, could be applied to dark soliton prop- 
agation |lf)j | or low-dimensional quantum degeneracy in- 
cluding the Tonks gas regime Mil. H^ . Toroidal traps 
have been around for some time, being employed in early 
experiments with Sodium vapors |l3l | . In this case a blue- 
detuned laser was used as a measure to counteract Ma- 
jorana spin flips; a loss mechanism which can be prob- 
lematic in the cooling required for condensation. In more 
recent experiments the toroidal topology has been used 
in the study of vortex nucleation and superfluidity [3| ■ 

Recent theoretical work concerning toroidal conden- 
sates has concentrated on excitations in traps which 
have some time dependence in their topology |15| and 
on vortex-vortex dynamics which have been shown to be 
strongly distorted by such a geometry 0] . The stability 
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of multiply-quantized toroidal currents has been studied 
by |{| [l5| . The spectrum of single-particle excitations for 
cigar-shaped toroidal traps with circulation has also been 
considered ,4f . Although a wide variety of toroidal trap 
parameters is possible, one of the advantages of such a 
system, the most easily accessible experiments appear to 
be based on oblate shaped traps. In this paper we study 
results for the spectrum of collective excitations of oblate 
toroidal condensates within the Bogoliubov approxima- 
tion, and explores the dynamical stability of ring currents 
around the torus. The main features we note are generic 
to this design of trap and would apply to similar geome- 
tries. Perturbations of superfluid ring currents by an off- 
set of the trapping potential are studied. For example, in 
a toroidal trap the central repulsive potential acts a pin- 
ning site for the vortex and thus the stability of the flow 
can be studied. The excitation spectrum, mode densi- 
ties, flow rate and center-of-mass motion for this system 
are obtained by employin g bo th a time-dependent and 
time-independent method |l7l l!s[. A simple, but accu- 
rate, formula is presented for the lowest angular acoustic 
modes of excitation, and the splitting energy of these 
modes when a background current is present. 



II. THE MODEL 

A. The weakly-interacting Bose gas 

For a cold dilute weakly-interacting gas, the ground 
state (condensate mode) dominates the collective dynam- 
ics of the system. In experimental realizations one can 
achieve temperatures such that T <C T c (typically 0.1 
to 1/xK) and densities such that the gas is weakly in- 
teracting and highly dilute. Under such conditions, the 
condensate of No ^> 1 atoms is well described by a mean 
field, or wavefunction, governed by the Gross-Pitaevskii 
equation, and the quasiparticle excitations are acoustic 
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waves within this field. If the perturbations of the con- 
densate are small, then it is appropriate and convenient 
to use the linear response approximation, which is equiv- 
alent to the Bogoliubov approximation for single-particle 
excitations in highly-condensed quantized Bose gases at 
zero temperature. 

Consider a dilute system of Nq atoms, each of mass 
m a , trapped by an external potential V ex t(x, t) and inter- 
acting weakly through the two-body potential V(x,x'). 
At low temperatures and densities, the atom-atom inter- 
action can be represented perturbatively by the s-wave 
pseudopotential: V(x,x') = {Ai:h 2 a s /m a )5^ (x — x'), 
and a s is the s-wave scattering length. The dynamics 
follow from the Hartree variational principle: 



0.008 



dt 



d 3 a 



ip*[H + ±gip*ip -iHdtltp = (1) 



where g = (4ttH 2 /m a )N a s , H a = -(h 2 /2m a )V 2 + V cxt - 
H, and the chemical potential fj, plays the role of a La- 
grange multiplier. Supposing that the temperatures are 
sufficiently low that the condensate can be represented 
by a Bogoliubov mean field (f>. The single-particle exci- 
tations can be described by the linear response ansatz: 



ip{x,t) = a (t)(f>(x) 



(2) 



+ a*{t)v*{x)e 



-\-iujjt~\ 



where (f> represents the highly-occupied condensate; that 



is, |ao| ~ v N ^> \a,j\, j ^ 0. From the variation 5(f)*, and 
linear expansion in the small parameters a,j , a* taken as 
constant, the stationary Gross-Pitaevskii equation and 
Bogoliubov equations follow: 



H^ + g\4>\ 2 4> = {) 



(3) 



with (») = 1 and j d 3 x f(x)*g(x) = (f,g). The Bo- 
goliubov modes are solutions of the coupled linear equa- 
tions: 



(H + 2g\ ( f>\ 2 )u J + g<f 2 v 
{H a + 2g\cj ) \ 2 )v 3 +g^ 2 u 



+hujjUj 
-hw 3 v 3 



(4) 
(5) 



Time-reversal symmetry of equations (|4I5|) is reflected 
in the fact that every set of solutions {u>j , Uj , Vj } has 

, u*j} and the normalization 



a corresponding set {- 
convention is: {u,i,Uj) 



- (vi,Vj) = 6 i: 



B. Toroidal condensates 

Conventional atom traps provide a confinement of 
the condensate in the radial and axial directions. We 
write this as a spheroidal potential of the type: V — 
^m a u> 2 (r 2 + A 2 z 2 ), where r is the radial coordinate, u) r 
the radial frequency, and z the axial coordinate. Here 
A is the aspect ratio of the harmonic potential, so that 
A > 1 flattens the condensate into an oblate pancake 
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FIG. 1: Torus potential (dashed line) and condensate den- 
sity </>| 2 (full line) as a function of radius in the the plane 
z = 0. The radial potential corresponds to /io = 50hur, and 
a = 2^/ 'h/2rn a LOr. The condensate circulation is n — and 
the interaction strength is C = 1000. The scaled condensate 
probability density |</>| 2 is given in units (fi/2m n u r ) 3 ' 2 



shape. For convenience we use scaled dimensionless units 
for length, time and energy namely: (h/2m a uj r ) ? , uj^ 1 
and huj r , respectively. Typically, Lu r ~ 2tt x 200Hz, so 
that the natural timescale for oscillation is of the order 
of milliseconds. If this potential is supplemented by a re- 
pulsive core then the toroidal shape can be realized. The 
conventional method is to use blue-detuned laser light 
to create a dipole force. In the experiment of Davis et 
al. the AC Stark shift from the green light of an 
argon-ion laser (A = 514nm), detuned from the D-linc 
(A = 589nm) of the trapped sodium atoms created the 
dipole force. If the laser beam is aligned along the 2-axis 
of the trap at the diffraction limit focus then the poten- 
tial if a function of r only: U(r) = hl(r)(8r 2 A./ sot ) _1 , 
where I sa t is the saturation intensity of the level, A is 
the detuning and r the lifetime of the upper atomic level 
[lflj . Across the focus the intensity variation is Gaus- 
sian: I(r) = Iq exp(— r 2 /a 2 ), where the on-axis inten- 
sity is related to the power P and waist of the beam 
wq by Iq = 2P/ttwq. So the width and height of the 
central barrier can be controlled by the these parame- 
ters. Consider a typical case for a trap containing sodium 
atoms; an argon-ion laser beam of 3.5 W focused to a 
waist wo = 30/im, would give a frequency shift of 7 
MHz at r — 0. Compared with trap frequencies of order 
uJr/2-K ~ 100 — 400Hz, this would be sufficient to create a 
hole along the trap axis and form a toroidal condensate. 
The combination of the fields gives the external poten- 
tial: V; xt (x,0) = \m a u 2 r {r 2 + \ 2 z 2 ) + h exp(-r 2 /a 2 ), 
where the repulsive central barrier is defined by height 
and width parameters h$ and a, respectively. The radial 
minimum of the well is displaced to tq — (T^/ln(4/io/c 2 )- 
For convenience we use scaled dimensionless units for 
length, time and energy namely: (fi/2m„w r ) 2 , lu^ 1 and 
ftto ri respectively Then it is convenient to define the in- 
teraction strength bv the dimensionless parameter C = 
8TrNoa s (h/2m a uj r )~2, so that as millions of atoms oc- 
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cupy the condensate No — ► oo,C — > oo, while the ideal 
gas corresponds to C — > 0. For the pancake geometry 
A> 1 the central barrier height ho controls the transition 
from spheroid to torus. The radius ro of the potential well 
minimum provides a guide to the size of the torus. Col- 
lective excitations for a quartic toroidal potential of the 
form: V(r, z) ex t = \ m a^Mr ~ r ) 2 + X 2 z 2 } were consid- 
ered by Salasnich et al Q and the variation in chemical 
potential and frequency with respect to particle number 
was studied. The geometry chosen £| was a prolate shape 
such that A = 1/ v8 in contrast to the oblate case consid- 
ered here. An example of a toroidal condensate density 
in static equilibrium is shown in figure ^ The trap pa- 
rameters are ho — 50hio r , and a = 1\J h/2m a uj r so that 
ro ph 4, and the interaction strength is C = 1000. For 
strong interaction, high C, the radius of gyration (tk) 
about the symmetry axis is a better estimate of conden- 
sate radius; where r\ = J \4>\ 2 r 2 d 3 x/ J \<fi\ 2 d 3 x. For the 
pancake shape (ho = 0) rx = 3.2, but as the barrier is 
raised to ho = 10 then there is some depletion of conden- 
sate at the center, so that it expands to tk = 3.8. Finally 
for a high central barrier, ho — 50, the condensate is ex- 
cluded from the trap axis and a narrow ring is formed 
with Tk = 4.6 as shown in figure ^ Here the maximum 
condensate density is located close to the potential energy 
minimum ro w 4 in accord with hydrostatic equilibrium 
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FIG. 2: Excitation spectrum for the axisymmetric modes 
of semi-toroidal condensate without circulation (k = 0) as 
a function of interaction strength C. Trap parameters are 
A = y/7, (7 — 2. The symmetries of the modes are: m = 
(o) ,\m\ = 1 (o) , |m| = 2 (*) . 
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FIG. 3: Excitation spectrum for the axisymmetric modes of 
toroidal condensate without circulation (k = 0) as a function 
of interaction strength C. Trap parameters are A = V7, a = 2 
and ho = 50. The mode symmetries are: m = (o),|m| = 
1 (o), \m\ = 2 (*) 



will be presented for the first few monopolc, dipole and 
quadrupole excitations only. For finite C, the spectrum 
of excitations must be determined by numerical solution 
of equations (|3I4|I and ©. Separating variables gives: 



(r,z,<p) = 4> K (r,z)e iK * 



(6) 



so that the condensate, with circulation k and real am- 
plitude 4>, is the solution of the equation: 



(7) 



where 



£( S ) = — 



2m„ 



d 2 

dr 2 



1 d d 2 

r dr dz 2 



\m a u? r (r 2 + \ 2 z 2 ) + h exp(-r 2 /(T 2 ) + 2g\4> K \ 2 - f i 

The cylindrical symmetry of the condensate means that 
small amplitude excitations can be described by radial 
(n r ), axial (n z ) and rotational (m) quantum numbers, 
each with an associated parity. Here m — ng — k = 
0, ±1, ±2, . . . denotes angular momentum with respect to 
the condensate, while n r , n z — 0, 1, 2, 3, ... . The quasi- 
particle amplitudes: 



Vn r ,m 



t ^(r,z,ip) = u nr , nz (r,z)e^ m+K ^ (9) 
^(r,z, V ) = v rir ^(r,z)e^^ (10) 



C. Bogoliubov spectrum 

In this paper we are primarily interested in the low- 
lying axisymmetric excitations arising from weak per- 
turbation of the condensate, so that detailed results 



with corresponding angular frequency, u) ri 
lutions of the eigenvalue problem 



,„ 2 , are so- 



£{m+K)u{r, z) +g 4>l v(r, z) = +huu(r,z) (11) 
C{ m - K )v(r, z) +g 4>l u(r, z) = -huv(r,z) (12) 



The equations are discretised on a 2D-grid using La- 
grange functions 0; the radial coordinate is defined at 
M grid points (n, r-2, . . . , Tm) and the axial coordinate 
at N points (zi, 22, ... , 2/v). Therefore: 

M N 

k(r,z) = ^2J2^^ z X k l/2 ^ 1/2 h(r)m(z) (13) 

k=l 1=1 

where /, g are Lagrangian interpolating functions such 
that 



f*( r )fk(r) 2irr dr w A 4 (5 ife 



'■(z)gi(z) dz w Aij<Sj-; 



(14) 
(15) 



The functions for the r-coordinate are chosen to be gener- 
alized Laguerre polynomials , scaled to encompass the 
entire condensate, with typically M — 50 mesh points. 
On the other hand, Hermite polynomials are used in the 
z-direction so that 



N-l 



(16) 
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FIG. 4: Excitation spectrum for modes of a condensate 
with circulation k = 1. The mode symmetries arc \m\ — 
0(o), l(o), 2(*). Frequencies are plotted as a function of 
barrier height ho in a trap characterized by A = \fl and 
a = 1\J h/2m a uJr, with interaction strength C = 1000. See 
also table HTl 



where xi( z ) = h N 2 w(z)^ Hi(z). and Hi(z) are the Her- 
mite polynomials associated with weights w(z) = er^ 
and normalization factor /ijy — 2 N n 1 / 2 N\. A high de- 
gree of accuracy was found with only N = 30 points. 
The resultant eigenvalue problem was solved using the 
following approach. The condensate density was found 
using Newton's method for equation Q. The decoupled 
linear eigenvalue problem I|11I12(I was solved by conver- 
sion to Hessenberg QR form followed by inverse iteration 
plj . Convergence was established by a combination of 
grid scaling and number of mesh points, so that at least 
six-figure precision was assured for all frequencies (see 
table |H and table HI. 



D. Spectrum of torus excitations 

The Bogoliubov mode frequencies for a toroidal con- 
densate without circulation (k = 0) arc shown in figures 
121 and |3| for two different barrier heights. The varia- 
tion in frequencies as the particle number (C) increases 
is shown. In the absence of circulation for the low barrier 
ho = 10 (figure I2J) the effect of increasing the interaction 
strength, or equivalently increasing the number of atoms 
in the condensate, is to spread the spectral lines of the 
low-lying modes, while the highly-excited state frequen- 
cies are relatively insensitive as C increases. In contrast, 
to the prolate toroid results Q| the effect of increasing 
atom number, increasing C, leads to an increase in the 
gap between the low-frequency excitations. The three 
to = modes shown are, in increasing frequency, the 
gapless mode u> = 0, the first breathing mode, and the 



TABLE I: Angular frequencies u/ui r of the first-excited ra- 
dial dipole modes. The table compares the numerical accu- 
racy of the time-dependent linear response method (TDLR) 
with e — 0.02 with the time-independent methods (BdG) and 
agreement is better than 2% in all cases. The condensate in- 
teraction strength is C = 8TrNoa 3 (h/2rn a uj r )~ 4 = 1000 and 
the axial mode is n z — 0. Results are presented for k = and 

K=l. 



mode 


ho / hu> r 
TDLR 


= 10 
BdG 


ho /huj r 
TDLR 


= 50 
BdG 


K = 
m = +1 


0.772 


0.771 


0.563 


0.554 


K = 
m — —1 


0.772 


0.771 


0.563 


0.554 


K = 1 
m = +1 


0.875 


0.874 


0.657 


0.652 


K = 1 
m — —1 


0.496 


0.501 


0.448 


0.447 



first axial dipole mode uo — \fl . The lowest degenerate 
modes m = ±1 are the fundamental radial dipoles, at 
higher frequency lo ~ 2 we observe an octupolc mode, 
and finally the xz-quadrupole at lo ~ 2.8 is very close 
to the axial dipole. The increasing density of states for 
higher excitations is familiar from studies of spherical 
condensates and is reproduced here. Finally the lowest 
pair of to = 2 states, corresponding to quadrupole ex- 
citation in the radial coordinate are far below the next 
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TABLE II: Excitation frequencies uj/uj r of the toroidal con- 
densate with circulation k = 1 with variation in the cen- 
tral barrier height ho- The condensate interaction strength is 

C = 8nN as{h/2m a LJr)~^ = 1000 



m 


ho/fUOr = 


ho/hulr = 10 


ho/haj r = 50 






9 f\Af\ 






3.310 


3.243 


3.954 


m = 


3.471 


3.600 


3.970 




4.691 


4.480 


4.771 




4.728 


4.807 


5.609 




4.936 


5.027 


5.882 








n 447 


771 = 1 


1.000 


0.875 


0.652 




2.613 


2.417 


2.451 




3.085 


2.942 


2.823 




1.197 


1.061 


0.890 




1.640 


1.498 


1.260 


\m\ = 2 


3.499 


3.037 


2.790 




3.322 


3.321 


3.114 




4.925 


4.530 


4.396 




4.840 


4.561 


4.471 



higher excitations of this symmetry. 

As the toroidal shape is more sharply defined (ho = 50) 
the high-frequency modes become more tightly grouped 
(figure |2J) and the gap to the lowest-order excitations 
widens. The long-wavelength low-energy modes are cir- 
culations around the torus, as shown in figure [3] The 
acoustic modes 20] appropriate in the limit C — > oo, are 
sinusoidal oscillations in density 8p and follow the micro- 
scopic wave equation, in our units: 

-u; 2 5p = 2V • ([fx - V cxt ]VSp) (17) 

In the Thomas-Fermi limit, the condensate occupies the 
region Q where fi — V cx t > 0. The angular waves satisfy 
the periodic boundary conditions 5p((p + 2i:m) — 6p((p). 
Taking the volume average of equation (|17fl over the ra- 
dial and axial density gives: uj — \m\uj s where 

cu 2 I 2-Krdrdz = 2C r - 2 (j> 2 2TTrdrdz (18) 
Jn Jn 

So that: 

u = \m\y/2C(4?r- 2 ) (19) 

The results in figure [2] reflect the change in dispersion 
relation u> ~ m 2 to to ~ m for the lower frequency modes 
as the interaction strength increases. The value of lo s 
varies slowly with changes in the interaction strength C, 
see figure |21 It is worth noting that the simple formula 
(|17|l gives a good estimate of the frequency spectrum. 
For example, for C — 1000 u> s w 0.57 in good agreement 
with the results shown (figure while for C = 4300 the 



value is ui s w 0.78. The relation (|19(l is quite general 
and can by applied to any toroidal geometry. For a large 
radius but tightly confined torus, so that (r~ 2 ) w r^ 2 
the relation simplifies to the ID result as expected: 

u = \my2C/VL r^l (20) 

where the wavenumber is quantized according to \m\/rK, 
and the speed of sound is proportional to the square root 
of the atom density. One might assume, for a given radial 
and axial mode, that this pattern of frequencies for the 
m sublevels would be repeated for the higher frequency 
bands. In fact this is not the case as confirmed in figures[21 
and|3| The bands become mixed and the u> ~ |m| relation 
does not hold. A more extensive list of the frequencies 
including the highly excited modes is given in Table ITT1 
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FIG. 5: Center-of-mass trajectory of a toroidal condensate 
current k = 1 displaced from the trap center. The torus has 
parameters A = v7, o — lyj h/2rn a LO r and ho = 50hu> r and 
interaction strength C = 1000. The figure shows the path 
traced out over a total time T = 100 in the horizontal plane 
by the condensate center-of-mass following a small horizontal 
displacement of the trap axis e = 1. The combination of the 
x-dipole and y-dipole oscillations leading to precession of the 
condensate. These dipoles modes are degenerate and nearly 
in phase. The spectrum density of the a;-dipole oscillations is 
shown in figure |7| 



E. Toroidal current flow 

If the condensate itself has a current c / 0, then the 
central barrier acts as a pinning site for the vortex. Fur- 
thermore the ±|m| degeneracy of the excitation modes is 
removed. Perturbations of the external potential are of 
interest since it allows us to study imperfections in the 
toroidal trap potential, and also the dissipation of the 
ring current. We propose that these modes can be stud- 
ied by a lateral displacement of the trap. In Figure 
the spectrum of excitations of a singly-quantized current 
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FIG. 6: Center-of-mass trajectory of a toroidal condensate 
(current k = 1) after radial displacement. As in figure |S] 
the torus has parameters A = VT, o — 2sJ h/2m a Lo r and 
ho = 50huu r , however the atom number is much higher in this 
case:C = 4300. The figure shows the path traced out over 
a time T = 100 in the horizontal plane by the condensate 
center-of-mass following a small horizontal displacement of 
the trap e — 1. The center of mass performs a circular motion 
due to two effects: the x and y dipole states are not degenerate 
in frequency and the amplitudes are out of phase, by roughly 

~7f/2. 



loop k = 1 are presented. In this figure we consider the 
transition from a conventional pancake geometry ho = 0, 
through a semi-toroidal trap ho = 10, before reaching 
the torus shape ho = 50. The presence of the barrier 
creates a narrower radial well, so that the m — ra- 
dial excitations are pushed to higher frequencies. The 
axial m = dipole mode u> — A w 2.65, is unaffected of 
course. The splitting of the m-degeneracy is greatest for 
the spheroidal ho — case, in particular the lowest pair 
of \m\ = 1 modes. The co-rotating mode m = +1 has a 
positive frequency lying close to the lowest m — 2 mode. 
The anomalous or counter-rotating mode m = — 1 has a 
negative frequency, but positive norm. Anomalous modes 
arise when a condensate has a stable topological excita- 
tion such as a vortex, indeed the presence of a vortex is 
a necessary condition for an anomalous excitation as its 
presence is required to satisfy conservation of energy and 
momentum in the system [23. In this case the central 
barrier raises the anomalous mode to positive frequen- 
cies effectively creating a stable pinning site for toroidal 
flow. The stability of high circulation current in toroidal 
traps has been discussed in detail by Busch et al. [l5| . 
The splitting of the \m\ = 1 pair reduces as the barrier 
rises and the condensate becomes toroidal. The closing 
of the energy gap reflects the fact that the barrier expels 
condensate density from the center so that the m = — 1 
cannot occupy the vortex core region and create a large 
energy gap. The |m| = 1 modes will be evident when 
dipole excitations of the circulating flow are discussed. 



The m = ±2 states, degenerate in figure|3], are also split 
by the condensate flow and in the toroidal limit, ho = 50. 
The energy splittings are proportional to \m\ the current 
momentum. This follows from the Bogoliubov equations 
when the centrifugal energy terms can be considered as 
perturbations. First-order perturbation theory applied 
to equations <|ll|l and (|12f> gives the energy shift: 



h 2 (m + n) 2 _ 2 

Ahw w (u ,r u ) 

2m„ 



h 2 (m — k) 
2m a 



■(v ,r 2 v ) 
(21) 

where uo and vq are the unperturbed quasiparticle states. 
This gives the energy splitting AE — h(iv + — cj_): 



AE w i\m\n x 



2m a 



(mqU - VoV )r 2 d 3 i 



(22) 



If the condensate and excitations are tightly confined 
around the radius tk then: 



AE w 4|to|k x 



2m a r 2 K 



(23) 



The splitting decreases as C, and hence tk increases. 
Referring to TableHlfo 1 ' ho = 50, the |m| = 1 splitting is 
close to that given by the perturbation expression (equa- 
tion l23fl AE w 0.19. However the lowest \m\ = 2 modes 
in Table[l]have a splitting (0.37) very slighty smaller than 
the pertubation result 0.38. 

The simplest experimental scheme to observe these 
modes is to displace the potential in the horizontal plane 
and follow the motion of the condensate center-of-mass. 
Consider a sudden adjustment of the trapping potential 



V oxt (x, y, z) -> Fcxt (x - e, y, z) 



(24) 



The effect on the center-of-mass is shown in figures |5] 
and figure[S]in which the trajectory is traced out as time 
evolves. These graphs were obtained by direct solution of 
the time-dependent Gross-Pitaevskii equation, following 
from taking arbitrary variation of ip* in l|T)l: 



[Ho + 9 \iA 2 - ihd t ] v> = o 



(25) 



For very small perturbations the nonlinearity is small and 
the time-dependent linear response method (TDLR) is 
equivalent to the stationary Bogoliubov equation method 
outlined above. This employs a direct numerical solution, 
in combination with spectral analysis can also be used to 
determine the frequency spectrum 01 an d the strength 
of each frequency component, that is the population of 
modes and density fluctuations of the collective excita- 
tions. This can be done efficiently and accurately using 
spectral methods. The initial state can be determined 
by imposing the phase according to the value of k and 
evolving in imaginary time The calculations pre- 

sented here were performed on a 64 x 64 x 64 grid. The 
scaling of grid spacing according to the degree of confine- 
ment ensures maximum spatial resolution of the conden- 
sate, so that the z dimension is more tightly confined. 
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FIG. 7: Spectral density of s-dipole of a condensate with 
circulation k = 1 corresponding to figure Px (u) The torus 
parameters are A = \/7, c = 2^/h/2m a uj r and fao — 50&j r . 
The condensate interaction strength is C — 1000 with trap 
displacement e = 1. The figure shows the dominance of 
three modes, the low frequency m = ±1 doublet, and an 
excited m = — 1 mode. The center of the doublet corre- 
sponds to the frequency u>s = ^/2C(j — 2 ) ~ 0.58. In the 
graph the frequency peaks are shifted from the linear re- 
sponse results given in Table U] and have a broader splitting 
than the Bogoliubov approximation. The a;-dipole modes, 
shown above, have frequencies u) x = 0.421, 0.666, 2.237, 2.441 
while the corresponding y-dipole modes have frequencies: 




FIG. 8: Spectral density of :r-dipole of condensate with cir- 
culation k = 1 corresponding to figure |S| P m (ui) The torus 
parameters are A = \/7, a = 2^/h/2m a uj r and ho = 50hco r . 
The condensate interaction strength is C = 4300 with trap 
displacement e = 1. The center of the doublet corresponds 
to the frequency cu s = ^2C{cj> 2 i — 2 ) ~ 0.78 The figure shows 
the shifting upwards of the mode frequencies and a narrow- 
ing of the m = ±1 doublet splitting consistent with an ex- 
panded condensate with a higher sound speed. The modes 
shown above have frequencies: ui x = 0.625, 0.804, 2.095, 2.284. 
The corresponding y dipole modes, not shown, have similar 
strengths but with frequencies u y = 0.657, 0.792, 2.086, 2.275. 



The time step chosen, At w 0.001 x 2ir/oj r , is dictated by 
the characteristic timescale of the excitations. A power- 
spectral-density estimate which used a Kaiser windowing 
function was employed to analyse the various moments of 
interest in terms of its component frequencies. For exam- 
ple the x-dipole moment can be transformed as follows: 



e- iult (x)dt 



(26) 



The spectral data obtained for a very small trap dis- 
placement e = 0.02 confirms the mode frequencies from 
the Bogoliubov-dc Gennes equation, and establishes the 
accuracy of the method. As shown in tablc[IJ the frequen- 
cies using this method are accurate to within a few per 
cent. However a realistic measurement process based on 
density imaging requires a much larger amplitude mo- 
tion in order to resolve the fluctuations. In this case, 
the nonlinear effect can be important. In our scheme, 
which we propose as a viable method to measure the 
spectrum, we consider the trap displaced by a substan- 
tial amount e = 1 in order to resolve the oscillations in 
the plane. When circulation of the toroidal current is dis- 
turbed in this manner the effect is to create large counter- 
circulating currents. If these currents are in phase and of 
equal magnitude the condensate would execute a pedu- 
lum motion. In general the motion will be 2D if the 
symmetry is broken. In terms of experimental observ- 
ables, the effect appears as oscillations of the center of 
mass with a precessional motion. Similar precession mo- 
tion for quadrupole excitation of a spheroidal vortex state 
has been noted ^3 • The locus of the center-of-mass in 
the horizontal plane from t = to t = 100, is shown in 
figure[S]for C — 1000 and ho/hujo = 10. Initially the con- 
densate moves from x = Q, y = towards the equilibrium 
point x = 1,2/ = 0, however this motion is converted an 
irregular precessing pendulum motion. The rate of pre- 
cession is proportional to the splitting of levels. However 
the center-of-mass motion is not simple. A clearer under- 
standing of the motion follows from a spectral analysis of 
the x and y dipole moments. In fact, the apparently irreg- 
ular motion is dominated by only a few components the 
low frequency dipoles. In both the case ho = 10 (figure 
0) and ho = 50 (figure |SJ the m = +1 mode is stronger 
than the mode m = — 1 and this dictates a motion with 
positive helicity, that is the condensate precession is in 
the same sense as the flow. However the frequencies are 
significantly different from the linear model and the x- 
y degeneracy is removed due to the nonlinearity of the 
response in the x direction. In figure the torus pa- 
rameters are A = v7, & = 2y/h/2m a u) r and ho — 50fku r . 
The condensate interaction strength is C = 1000 with 
trap displacement e = 1. Comparing the x frequencies 
with the values of table [I] we note the y-dipole frequen- 
cies are very similar in both cases, but the x-dipole pair 
splitting is increased when the trap displacement is this 
large. The motion shown in figure [S] is a superposition of 
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a faster x dipole with a slower y dipole. The results for 
the higher interaction strength C = 4300 show the split- 
ting of the |m| = 1 lines narrowing as the size of the con- 
densate increases. The figure shows the shifting upwards 
of the mode frequencies and a narrowing of the to = ±1 
doublet splitting consistent with an expanded conden- 
sate with a higher sound speed. The modes shown have 
frequencies: w x = 0.625, 0.804, 2.095, 2.284. The corre- 
sponding y dipole modes have similar strengths but with 
frequencies u y = 0.657, 0.792, 2.086, 2.275. The nonlin- 
ear interactions give different splittings for the x and y 
dipoles: 5u) x = 0.179 compared with 5u> y — 0.135 both 
of which are slightly smaller than the linear perturbation 
theory (|23|l 6oj w 0.200. The major contrast between 
figures and El is the orbiting motion observed for the 
higher value of C, although the frequency components 
shown in figures and [5] are of similar relative strengths, 
the x and y dipoles have a large phase difference close to 
7r/2 for C = 4300 and explains the more circular char- 
acter of the motion. This is slightly surprising in that 
the usual behavior of collective mode of oscillations is 
that the frequency dependence on TV is relatively weak 
for large N. In this respect, collective excitation frequen- 
cies are usually a poor method of estimating condensate 
fraction compared with the hydrostatic pressure. How- 
ever, in this case the toroidal geometry seems to retain 
the quantum features of the excitations even for large 
particle number. 

III. CONCLUSIONS 

In this paper we studied the spectrum of collective exci- 
tations of oblate toroidal condensates within the Bogoli- 



ubov approximation, and explored the dynamical stabil- 
ity of ring currents around the torus. The main features 
we noted are generic to this design of trap and would ap- 
ply to similar geometries and would be exhibited in ex- 
periments on toroidal condensates. The transition from 
spheroidal to toroidal geometry of the trap displaces the 
energy levels into narrow bands. The lowest-order modes 
are quasi one-dimensional circulations with dispersion re- 
lation uj ~ muj s with to quantized. When the condensate 
has an toroidal current flow, the low-energy circulations 
are split into co-rotating and counter-rotating pairs. A 
simple, but accurate, formula is presented for the lowest 
angular acoustic modes of excitation, and the splitting 
energy when a background current is present. The sta- 
bility of the flow is maintained by the central barrier 
acting as a pinning site. 



Instabilities will most readily occur by radial dipole ex- 
citations. The dipole states are nearly degenerate in the 
toroidal limit. Small condensate currents when disturbed 
create slow precessional motion through the center of the 
trap through mixing of the degenerate x and y dipole 
states. For a tightly confined torus, these modes are well 
separated in energy and accessible to observation. We 
propose an experiment that would detect these modes. 
A realistic measurement requires resolution of the global 
condensate motion. Large displacements, such as those 
we propose, give rise to small nonlinear shifts in the fre- 
quencies beyond the Bogoliubov model. We also note 
that the highly occupied condensate performs a circular 
orbiting motion of the center-of-mass. Such features are 
resolvable and experimentally measurable with current 
technology and expertise. 
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